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Abstract The standard text-book Jacobi equation (equation of geodesic 
deviation) arises by linearizing the geodesic equation around some chosen 
geodesic, where the linearization is done with respect to the coordinates and 
the velocities. The generalized Jacobi equation, introduced by Hodgkinson in 
1972 and further developed by Mashhoon and others, arises if the lineariza- 
tion is done only with respect to the coordinates, but not with respect to 
the velocities. The resulting equation has been studied by several authors in 
some detail for timelike geodesies in a Lorentzian manifold. Here we begin by 
^ I briefly considering the generalized Jacobi equation on afhne manifolds, with- 

out a metric; then we specify to lightlike geodesies in a Lorentzian manifold. 
We illustrate the latter case by considering particular lightlike geodesies (a) 
, in Schwarzschild spacetime and (b) in a plane-wave spacetime. 

Keywords general relativity • light rays • Jacobi equation 



1 Introduction 



■ The Jacobi equation, also known as the equation of geodesic deviation, de- 

, scribes geodesies in the neighborhood of a reference geodesic. More precisely, 

' the Jacobi equation is the linearization of the geodesic equation around the 

reference geodesic, where the linearization is done with respect to the co- 
ordinates and the velocities. Thus, the Jacobi equation describes geodesies 
which are close to the reference geodesic and whose velocities are close to the 
velocity of the reference geodesic. The Jacobi equation has important ap- 
plications to General Relativity: For timelike geodesies, the Jacobi equation 
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describes the relative acceleration of freely falling particles around a freely 
falling reference particle. This relative acceleration, which is determined by 
the curvature tensor along the worldline of the reference particle, can be in- 
terpreted as the tidal force produced by the gravitational field. For lightlike 
geodesies, the Jacobi equation determines the shape of light bundles around 
a reference light ray. Again, the expansion and distortion of the bundle is 
determined by the curvature tensor along the reference light ray. Detailed 
discussions of the Jacobi equation (or geodesic deviation equation) can be 
found in almost any text-book on General Relativity, see e.g. Synge [T9] . 
Misner, Thorne and Wheeler [14] or Hawking and Ellis [8j. 

If a geodesic is close to a reference geodesic but its velocity is not, the 
Jacobi equation does not give a valid approximation for it. Such a geodesic is 
properly described by an equation which is linearized only with respect to the 
coordinates but not with respect to the velocities. Such a generalized Jacobi 
equation was brought forward, for timelike geodesies in a general-relativistic 
spacetime, by Hodgkinson [3] . This equation describes tidal forces on particles 
whose relative velocity is not small. Applications to astrophysics were given 
by Mashhoon [T2lll3j. The generalized Jacobi equation was independently 
rediscovered by Ciufolini [6] . More recently, its implications were studied in 
a series of papers by Chicone and Mashhoon [31I111S]. 

In geometric terms, the ordinary Jacobi equation approximates the geo- 
desic flow on a tubular neighborhood of the reference geodesic in phase space 
(i.e., in the cotangent bundle over the base manifold). By contrast, the gen- 
eralized Jacobi equation approximates the geodesic flow on a neighborhood 
in phase space that is unrestricted in the fiber dimension. 

All the references quoted above consider the generalized Jacobi equation 
for timelike geodesies, which is physically of particular relevance because of 
its relation to tidal forces. In this article we want to discuss the generalized 
Jacobi equation for lightlike geodesies. Then it describes light rays around 
a reference light ray. The domain where the generalized Jacobi equation is 
applicable but the Jacobi equation is not encompasses all neighboring light 
rays that are close to the reference light ray but whose velocities are not. 

Although we are mainly interested in lightlike geodesies, it is worthwile 
to note that the geodesic equation and, hence, the Jacobi equation and the 
generalized Jacobi equation can be formulated on an affine manifold, without 
a metric. All that is needed is a connection. We begin by deriving the affine 
generalized Jacobi equation in Section [51 In Section [31 we specify to lightlike 
geodesies of a Lorentzian metric. In this case, the generalized Jacobi equation 
can be used to describe light bundles with arbitrarily large opening angle 
around a reference light ray. Of course, if the opening angle is large, the light 
rays will leave the neighborhood of the reference light ray soon, unless they 
are being refocused towards the reference light ray. So the generalized Jacobi 
equation can be used, in general, as a valid approximation only close to the 
vertex of the bundle. The general results of this section are illustrated by 
considering a bundle around a circular light ray at r = 3to in Schwarzschild 
spacetime. In Section[4lwe consider lightlike Fermi normal coordinates which 
allow to write the generalized Jacobi equation in terms of the curvature tensor 



3 



along the reference light ray. As an example, we consider a bundle of light 
rays in a plane-wave spacetime. 

Throughout we work with the generalized Jacobi equation in the sense of 
Hodgkinson, Mashhoon etc. This equation is linearized with respect to the 
coordinates but retains full dependence of the velocities. As an alternative, 
one may approximate with respect to the coordinates and the velocities up 
to common order N. This approach, which dates back to Bazahski [Tj, can 
be used to set up an iterative scheme. One begins by solving the ordinary 
Jacobi equation, then calculates the corrections due to second-order terms 
and so on. This method was applied to particle motion in Schwarzschild and 
Kerr spacetimes by Kerner et al. p^Ff]. 

This paper is dedicated to Bahram Mashhoon on occasion of his 60th 
birthday. I have very much profited from discussions with him and from 
reading his papers. In particular, I wish to thank him for introducing me to 
the subject of the generalized Jacobi equation. 



2 The generalized Jacobi equation in an afRne manifold 

An affine manifold is a manifold with an afline connection V. In local coor- 
dinates, the connection is determined by its connection coefficients -T/.'^, 

Va„9, = (1) 

Here and in the following, the dimension of the manifold is n, and we use the 
summation convention for greek indices running from 1 to rt. 

Whenever we have an affine manifold, we can consider the geodesic equa- 
tion 

The solution curves x^(s) of ^ are the geodesies (autoparallels) of the affine 
connection. As a geodesic remains a geodesic under affine reparametrization, 
s I— *■ as -|- 6, it is usual to refer to s as to an "affine parameter" . Note that 
only the symmetric part of Fj^^ enters into the geodesic equation. The anti- 
symmetric part 

T!ia - , (3) 

which is called the "torsion" , drops out from the geodesic equation. As in the 
following we are interested only in the geodesic equation, and in equations 
derived thereof, we may replace any connection by its symmetrized version. 
Therefore we will assume in the following that the connection is symmetric, 
rif^ = Fj^^, . (For the Jacobi equation and some of its generalizations in terms 
of a non-symmetric connection see, e.g., Swaminarayan and Safko ^iSj.) 
Now fix a geodesic X^(s), 

d^X^' , , , dX" dX'' 

We will call it the "reference geodesic" , and we will assume that it is known, 
i.e., that we have the X^^ explicitly as functions of the curve parameter s. 
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Then the geodesic equation for a neighboring curve x^{s) = X'^{s) + ^^(s) 
results from inserting x^{s) — X'^{s) + £,^{s) into ^ and subtracting ([4]), 



ds^ 



, f dX'^ dt" \ / dX" dt" \ 



""^^ ds ds 

With X'^(s) known, ([5]) is a system of second order ordinary differential equa- 
tions for the n functions ^''(s). It is the exact geodesic equation, expressed in 
terms of the coordinate difference £,^{s) with respect to the reference geodesic 
X^{s). 

If we linearize ([5]) with respect to but not with respect to d^f^ /ds, we 

get 

d^C'' u , . I dr dX-" dC dr \ 



a 



ds^ '"^ V ds ds ds ds 



This is the generalized Jacobi equation, on an affine manifold, in arbitrary 
coordinates. 

By contrast, if we linearize ^ both with respect to and with respect 
to d£^^ /ds, we get 

d^^f" „ , , dC dX" „ „ , . -r dX" dX" 

This is the ordinary Jacobi equation in arbitrary coordinates. To recover the 
standard text-book form, we have to introduce the covariant derivative along 
X{s), which is defined by 

for any rj" , and the curvature tensor 

R^,,, = d,r!/, - drn, + ri:^ ~ r!/^ r^, . (9) 

Then a straight-forward calculation reduces ([7]) to 

^^-0, (10) 

which is, indeed, the standard text-book form of the ordinary Jacobi equa- 
tion. 

The generalized Jacobi equation ([5]) is a second order non-linear ordinary 
differential equation for £,^{s). It is non- autonomous because the coefficients 
rif^{X) and drrif^{X) are functions of the curve parameter s. It gives a valid 
approximation for all those geodesies for which the ^^(s) are small, whereas 
the d£^^{s)/ds need not be small. Of course, if the d£^^{s)/ds are not small, the 
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S,^{s) will in general remain small only for a small interval of the parameter s. 
Such geodesies will leave the neighborhood of the reference geodesic soon and 
the generalized Jacobi equation will, in general, give a valid approximation 
for them only on a small interval of the parameter s. 

The generalized Jacobi equation may be modified in two ways. 

— One may prefer to choose a different, non-afhnc, parametrization for the 
(approximate) geodesies. 

— One may wish to choose special coordinates such that the generalized 
Jacobi equation takes a simpler form and that it can be easier compared 
for two different geodesies. 

Both is possible in an affine manifold, without refering to a metric. 

As to the first item, it is often convenient to use on a geodesic a (non- 
afhne) parameter u which coincides with one of the coordinates, say w = x". 
This is possible for all curves which are transverse to the hypersurfaces 
cc" = constant. Then the dependence of the coordinate on the parameter 
is known, x"'{u) — u, and the geodesic equation reduces to a system of differ- 
ential equations for the remaining coordinate functions (x^ (u), . . . , (u)) . 
The condition that two curves X{u) and X(u) + S,{u), both parametrized by 
u = x", are geometrically close is equivalent to the condition that the 
are small for i = l,...,(n — 1). By contrast, if an afhne parameter s is used 
it may be that not all of the £,^{s) are small even if the two geodesies are 
geometrically close. This happens if the affine parameter on one of the two 
geodesies lags behind or runs ahead of the other. 

Therefore we will now rewrite the generalized Jacobi equation with re- 
spect to the new curve parameter u — x"'. We first observe that the ^ — n 
component of the geodesic equation ^ can be written in the form 

and the remaining (n — f ) components as 

d fdudx^ \ , . dx" dx°' / du „ 

d;U*7)+^-^^)d^^(d;) (1') 

Here and in the following, latin indices take values 1 to (n — 1). 

Calculating the first term of (|12p with the product rule and inserting (fTTI) 
results in 

Together with the equation a;"(it) — u, determines the geodesies (or, 
more precisely, those geodesies that are transverse to the surfaces x" = 
constant) with the non-affine parametrization by u = x". (|13p is a system of 
second order ordinary differential equations for the (n — 1) functions x^{u). 
Clearly, these differential equations are non-autonomous because the coeffi- 
cients depend on u — x"'. Also, as a result of the reparametrization we have 
now got a system of equations that is cubic, rather than quadratic, in the 
velocities. 
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In analogy to our earlier procedure we search for solutions of (|13|) in the 
form x'^{u) — X (u) + (^f^ (u) , with a known geodesic, now reparametrized 
by M = x". Of course, as a;"(u) = X"-{u) = u, we have = 0. Inserting 

x^{u) = X'^{u) + into (UHl), and linearizing with respect to the C but 

not with respect to the d^'^/du, gives us the reparametrized version of the 
generalized Jacobi equation, 

d^C , , . ( dX" dr de df^x 
V du du du du J du 

_r" (X) + ^)(^ + ^)^ 

\ du du J \ du du J du 



-9.r" (x)e((^ + ^)(^ + ^)(^ + ^) ^ 0. 

\\du du J \ du du J \ du du J 

We now turn to the question of whether the generalized Jacobi equation 
can be simplified by choosing special coordinates. To that end we will use 
Fermi coordinates. Text-books on general relativity treat Fermi coordinates 
near a timelike reference geodesic (or, more generally, near a timelike refer- 
ence curve) in a Lorentzian manifold, see e.g. Synge |19j . However, it is fairly 
obvious that, by an analogous procedure, one can introduce Fermi coordi- 
nates near an arbitrary reference geodesic in an affine manifold, without a 
metric. The construction is as follows. 

Let an affinely parametrized geodesic X{s) be given. Choose an n-bein 
(i.e., n linearly independent vectors) at one point of the geodesic, such that 
the nth vector coincides with the tangent vector of the geodesic, and trans- 
port this n-bein parallelly along X{s). Denote the n resulting vector fields 
along the geodesic by Ei{s), . . . , En{s). As X{s) is an affinely parametrized 
geodesic, i?„(s) coincides with the tangent vector dX{s)/ds for all s. Let 
exp denote the exponential map of the given connection which, as indi- 
cated above, is assumed to be torsion-free. Then every point in a sufficiently 
small tubular neighborhood of the given geodesic can be written uniquely as 
ex^Pxiu) {^^Ei (u)). The n numbers {x^ , . ■ • , x^~^ , u) are the Fermi coordinates 
of this point. By definition of the exponential map, this means that the Fermi 
coordinates are determined by the following property: The point with Fermi 
coordinates (x^, . . . ,x'^~^,u) can be reached by following the geodesic with 
initial point X{u) and initial tangent vector x^Ei{u) up to the affine param- 
eter 1. By construction, the Fermi coordinate u coincides along the geodesic 
X{s) with the affine parameter s. Along neighboring geodesies, however, the 
Fermi coordinate u will give a non-affine parametrization in general. 

Note that Fermi coordinates are well-defined only on a (possibly small) 
tubular neighborhood of a chosen geodesic X{s). Farther away from X{s) 
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they need not exist, because there may be points that cannot be reached by 
a geodesic starting on X{s), and they need not be unique, because geodesies 
issuing from a point on X{s) may intersect. 

The crucial property of Fermi coordinates is that, in such coordinates, the 
connection coefficients r'^^ vanish on X{s). To prove this, we first observe 
that Va,,9^ = along X{s) because the Efj_ are parallelly transported, hence 
ri^^{X) = 0. As the connection is symmetric, this imphes that rj^^{X) = 0. 
What remains to be proven is that r^^{X) = 0. We use the fact that 

V(a,+a,)(a, + 9^) = Va^a, + Vo^d, + Vs^dj + Va^d., . (15) 

As for any [n — l)-tuple (a;^, . . . , x"^^) the integral curves of x^di that start 
on X{s) are geodesies, the left-hand side and the first two terms on the right- 
hand side of US]) vanish on X(s). Thus, US]) tells us that r^(X)-Frj^(X) = 0. 
As our connection is symmetric, this implies that, indeed, r^-{X) = 0. 

If written in Fermi coordinates, the generalized Jacobi equation (|14p sim- 
plifies considerably because of rif^{X) = 0. Moreover, the use of Fermi co- 
ordinates reduces the generalized Jacobi equation to a standard form which 
facilitates comparison of this equation for two different geodesies. In Sec- 
tion [4] below we will make this explicit for lightlike geodesies in a Lorentzian 
manifold, where the Fermi coordinates can be further specified. 



3 The generalized Jacobi equation for lightlike geodesies in 
arbitrary coordinates 

We now specify the results of the preceding section to the case that V is the 
Levi-Civita connection of a pseudo-Riemannian metric g — g ^^dx^ dx" ^ 

= ^ 5^" ( d,g,p + d.g.p ~ dpg^a ) , (16) 

where, as usual, 

9^" Qua ^5^,. (17) 

We assume that the metric has Lorentzian signature (-f ,,...,+, —), and 
we want to discuss the generalized Jacobi equation for the case of lightlike 
geodesies, 

dx^ dx" 
ds ds 

If the dimension n of the manifold is equal to 4, lightlike geodesies can be 
interpreted as light rays in a general-relarivistic spacetime. The mathematical 
results, however, hold for any n. 

Let us briefly recall the well-known fact that, if the connection coefficients 
rif^ are given in terms of a metric via (fT6|) , the geodesic equation ^ can be 
derived from a Hamiltonian. As a matter of fact, Hamilton's equations for 
the Hamiltonian 

H{x,p) ^ l-g^'''{x)p,,p,. (19) 
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take the form 

^=P., (20) 

^ = lid,9nix)PpPa. (21) 

The derivative with respect to s of (pij) . together with (PO)) and ([TB)l . indeed 
reproduces the geodesic equation By (PO)) . the side-condition takes 
the form 

i7(a;,p)=0. (22) 

A reparametrization of the solution curves can be achieved for lightlike 
geodesies in a particularly convenient way: We just have to multiply the 
Hamiltonian with an appropriate function of x and p. If we want to change 
from the afhne parameter s to the parameter u = x", we can proceed in the 
following way. We first recall that such a reparametrization is possible for 
all geodesies that are transverse to the hypersurfaces = constant. (This 
condition is true for all lightlike geodesies if and only if the hypersurfaces 
a;" = constant are spacelike. However, we may also consider the case that 
these hypersurfaces are lightlike or timelike; then the transversality condition 
only holds for some lightlike geodesies.) By (|20p and (|18p . the transversality 
condition is true if g"""^ {x)p^ ^ 0. This condition determines an open subset 
of the cotangent bundle (phase space). On this subset we can switch from 
the Hamiltonian to the modified Hamiltonian 

(23) 

Hamilton's equation with the Hamiltonian H read 

dx^ H{x,p)g^^^'{x) g^"'{x)p^ 



du g'^"{x)pa g''^"{x)pa ' 

dp^ _ Hix,p){dpg"P)ix)pp , {d^gn{x)ppPr 



(24) 
(25) 



du g"^'^{x)pcr 2g"''^{x)pa ' 

and the side-condition (|22p is equivalent to 

Hix,p)=0. (26) 
With the II — n component of ([M]) reduces to 



du 



1, (27) 



so the new curve parameter u coincides, indeed, with the coordinate x" (up 
to an additive constant that can be chosen at will). The other components 
of (gl) take the form 

dx' g^''{x)p^ 



du g^'^{x)pcr 



(28) 
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where, as before, latin indices range from 1 to (n — 1), and ([25|) reduces to 



Applying the derivative operator d/du to (pS)) . and using (^5]) . reproduces 
the equation (fTS]) for u-parametrized geodesies. We have thus proven that 
the lightlike geodesies of the metric g = g ^^{x)dx^ dx'^ , if reparametrized 
by M = x", can be derived from the Hamiltonian ([23| together with the 
side-condition (|^ . 

This gives us a convenient method of how to derive the generalized Jacobi 
equation for lightlike geodesies parametrized hy u — x^: Start out from the 
Hamiltonian ((23)) . Write Hamilton's equations with the side-condition ((26)) . 
Rewrite this as second-order equations for the a;^. Fix a solution X^^{u) to 
these equations, substitute x^{u) = +^^(it), and linearize with respect 

to the ^'^ but not with respect to the d^^^/du. 

The resulting set of equations describes lightlike geodesies X{u) + £^{u) 
which are close to the lightlike geodesic X{u) but whose tangent vectors need 
not be close to the tangent vector of X{u). For n = 4 these equations can 
be used, by choosing initial conditions appropriately, for describing a homo- 
centric light bundle around a central light ray X{u) in a general-relativistic 
spacetime. If the vertex of the bundle is at u = uo, say, we have to choose the 
initial conditions as Ci'^o) — 0. There are only (n — 2) = 2 independent solu- 
tions because the side-condition ((T8| fixes one of the d^"^ /du in terms of the 
others. This is in agreement with the intuitively obvious fact that a homo- 
centric light bundle has {n—2) — 2 dimensions transverse to the propagation 
direction. In contrast to the ordinary Jacobi equation, the generalized Jacobi 
equation can be used to describe homocentric bundles whose opening angle 
is arbitrarily large but the approximation is valid, in general, only for small 
values of w — uq, i-e., close to the vertex of the bundle. For larger values of 
u — uq, the will, in general, not be small, so the fact that we linearized 
with respect to these quantities may produce large errors. 

To be sure, there are special examples were the generalized Jacobi equa- 
tion holds for a large parameter interval. An example of this kind will be 
given in Sectional in this example even the Jacobi equation, the generalized 
Jacobi equation and the exact geodesic equation coincide. In general, how- 
ever, the generalized Jacobi equation for lightlike geodesies is a short-time 
equation, describing the temporal evolution of light bundles with arbitrarily 
large opening angles near their vertex. 

We illustrate the general results of this section with an example. 
Example 1: We want to calculate, with the help of the generalized Jacobi 
equation, the evolution of a light bundle around a circular geodesic at r = 3m 
in Schwarzschild spacetime. The Schwarzschild metric is 



du 



{d^,gP^){x)ppPr 
2.9"'^(x)p^ 



(29) 



p^jy dx^ dx' 



(1 



) df + 

r 



dr"^ 



+ dd^ + sin^i? dip^ . 



(30) 



2m 
r 



It is well-known that a lightlike geodesic that starts tangentially to the circle 
r = 3m, 1? = 7r/2, will stay on this circle. We want to write the generalized 
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Jacobi equation for lightlike geodesies near this circular geodesic. We will use 
the azimuthal coordinate ip for the parameter, u = = ip. This excludes 
geodesies tangent to a half-space (/? = constant. Therefore, our parametriza- 
tion allows us to treat bundles with any opening angle smaller than 7r/2 
around the circular geodesic, but not the limiting case that the opening an- 
gle is equal to 7r/2. Using for the parameter gives us directly the intersection 
of the bundle with any half-plane Lp = constant. 

After calculating the contravariant components g^'^ of the Schwarzschild 
metric, we write the Hamiltonian (|23p with x" = (/? , 

. 1 + + - ) . (31) 

Using the side-condition H{x,p) ~ 0, Hamilton's equations for the Hamilto- 
nian H take the form 

dip ^ dt r'svcL''dpt dt} sin^'dp^ 

du ^ du [r ~ 2m) ' du Pip ' 

dr r [r — 2m) sm^dpr 
du p^ 

dpip ^ p ^ = — = ^^^^Py C33\ 
du ^ du ' du smi} ' 

dpr {r — m) sink's p^ r^{r ~ 3m) sin^iJpj 

du p^p (r-2m)2p^ 

If we apply the derivative d/du to the expressions for dd/du and dr/du from 
([32]) . use (f33|) and the side-condition H{x,p) = 0, we arrive at the following 
second order system for and r{u). 

d^d 2cosd r d^\^ 



du^ sin I? 



/ dtf \ ^ 

— + sin-dcos-d. (34) 
\ du J 



d^r 2 cos -d d-d dr 2 / dr \^ , , //di9\2 „ 



du^- sin^ TuTu-^-KTu) +sin^.^j.(35) 
Obviously, this system of equations admits the solution 

r{u) 3m, d(u) = | . (36) 
In order to linearize around this circular lightlike geodesic, we write 

r[u) = 3 m -f i= r (") , m = I + e'C") ■ (37) 
V3 2 3m 

The numerical factors are chosen such that 

5(%,%) = = 1 (38) 
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dX r — 3m, d ~ 7r/2. Inserting ([37|) into ((34|) and ([35|) . and linearizing with 
respect to and ^""j but not with respect to d£f /du and d£^^ /du, gives us 
the generahzed Jacobi equation: 




du^ 9 m? 



(l)-l(l) « 



2^'' d^'' dc 2 ^ de 

9 du du 3-\/3' 



\ du ) 



If we hnearize also with respect to d^"^ /du and d^^/du, only the first term 
on the right-hand side of ([55)1 and of ([^0]) survives, 

This is the ordinary Jacobi equation. 

If we solve (j4T|) and (j42|) with initial conditions 

e{0) = 0, r(0) = 0, (43) 

d£'' dC 

-^(0) = -esinx, -^(0) = ecosx, (44) 

dip dip 

with X running from to 27r and e fixed, it gives us the cross-section of an 
initially circular bundle with opening angle proportional to e. Owing to the 
linearity of the ordinary Jacobi equation, the cross-section of such a bundle 
will be elliptic for all values of u — ip. As a, consequence of the plus sign on 
the right-hand side of (|42)) . in contrast to the minus sign on the right-hand 
side of (HU, the expansion of the bundle will increase in the direction 
and decrease in the direction; hence the major axis of the ellipse is in the 

direction. This reflects the fact that the circular geodesic at r = 3m is 
unstable with respect to perturbations in the direction but stable with 
respect to perturbations in the ^'^ direction. 

By contrast, if we solve the non- linear equations ((39|) and (f40|) with initial 
conditions and (|44p . the cross-section of the resulting bundle will not be 
elliptic. The larger the opening angle e, the stronger the deviation from the 
elliptic shape due to the non-linearities, see Figure [TJ 

This example demonstrates how the generalized Jacobi equation can be 
used for calculating the shapes of light bundles, beyond the small-angle ap- 
proximation that is inherent in the standard treatment based on the ordinary 
Jacobi equation. Of course, one has to keep in mind that the generalized Ja- 
cobi equation is, in genetral, a valid approximation only close to the vertex 
of the bundle. 
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Fig. 1 Cross section of initially circular light bundle around a circular lightlike 
geodesic in Schwarzschild spacetime, calculated with the generalized Jacobi equa- 
tion. The point {^^,^^) ~ (0,0) represents the circular geodesic at r = 3m and 
i9 = Tv/2. The axis is in the equatorial plane, pointing outwards. The axis is 
perpendicular 1;o the equatorial plane, pointing downwards. The picture shows the 
intersection with a half-plane ip = constant close to the vertex of the bundle, for 
four different opening angles. 



4 The generalized Jacobi equation for lightlike geodesies in Fermi 
coordinates 

We have already discussed how, on an affinc manifold, Fermi coordinates can 
be introduced near a reference geodesic X{s). The construction involved the 
choice of a parallely transported n-bein Ei{s),. . . ,En{s) along the chosen 
geodesic vifhich was arbitrary apart from the fact that En{s) should coincide 
with the tangent vector dX{s)/ds of the reference geodesic. If our connection 
is the Levi-Civita connection of a Lorentzian metric, we can further specify 
this n-bein. For a timelike reference geodesic, it is usual to require that the 
n-bein be orthonormal. Then Ei(s), . . . ,£"„_! (s) span the spacelike ortho- 
complement of dX{s)/ds at each value of s. Thus, in the Fermi coordinates 
{x^ , . . . the hypersurfaces u — constant intersect the timelike ref- 

erence geodesic X{s) orthogonally and are, therefore, spacelike near X(s). 
(Farther away from X{s) they need not be spacelike.) These are the Fermi 
normal coordinates treated in standard text-books on general relativity. In 
these coordinates the metric attains a standard form if written up to second 
order in the transverse coordinates. This standard form was derived by Man- 
asse and Misner [TT] and can be found, e.g., in Misner, Thorne and Wheeler 
[H, p.332. 
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Whereas this standard text-book treatment of Fermi normal coordinates 
assumes a timehke geodesic, here we are interrested in Fermi coordinates 
near a hghthke reference geodesic X{s). Then the n-bein cannot be chosen 
orthonormal because of the requirement that i?„(s) = dX(s)/ds. The best 
choice is to have £'„_i(s) hghthke with 

= -1, (45) 

and the remaining vectors i?i(s), . . . ,i?„_2(s) orthonormal and perpendic- 
ular to both En-i and En- (Condition (f45|) assures that En-i is future- 
pointing if En is future-pointing.) With this choice, the resulting Fermi co- 
ordinates (x^, . . . ,a;"~^,x"~^ — v,u) yield hypersurfaces u — constant that 
are lightlike where they meet the reference geodesic X{s). The hypersurface 
ti = is lightlike along the reference geodesic X{s) which is completely con- 
tained in this hypersurface. These lightlike Fermi normal coordinates were 
discussed in some detail in a recent article by Blau, Frank and Weiss [5]. 
In that article, the authors derive the general expression for the metric in 
lightlike Fermi normal coordinates, up to second order in the coordinates 
away from the reference geodesic. This result is the lightlike analogue of the 
above-mentioned Manasse-Misner representation of the metric near a time- 
like geodesic. If adapted to our sign and index conventions, it reads 

5^1^ dx^ dx'^ = —2dudv + Sab dx^ dx^ — Rinjn (u) x-' dv? (46) 

— - Rikjn (u) X-' dx^ du — — Rikji (u) x-' dx^ dx^ . . . 

Here the Fermi coordinates are denoted {x^, . . . , a;"~^, a;""^ — v,u), as out- 
lined above. As before, greek indices run from 1 to n and lower case latin 
indices k,l, . . . run from 1 to (n — 1). In addition, upper case latin indices 
A,B,... run from 1 to (n — 2). As usual, 6ab denotes the Kronecker delta. 
RpvaT{u) — gpf_i{u) R'^^^{u) is the purely covariant version of the curvature 
tensor Q, evaluated at {x^ = 0, . . . ,a::"~^ = 0,a;"~^ = v = 0,u), i.e., along 
the reference geodesic. The ellipses in (|46l) indicate terms of third and higher 
order with respect to the x*. 

With the contravariant metric components g'"' calculated, up to second 
order, from (|46l) . we can write the Hamiltonian ((23)) up to second order. This 
is enough to write the pertaining Hamilton equations up to first order, which 
together with the constraint (|26p will give us the generalized Jacobi equation 
in lightlike Fermi normal coordinates near an arbitrary lightlike reference 
geodesic. 

We will illustrate this with an example. 

Example 2: Consider a plane- wave spacetime in Brinkmann coordinates 
(x\ . . . ,x"''^,v,u), 

g^^dx^'dx" = -2dudv - hAsiu) x^ x^ du'^ + Sab dx^ dx^ , (47) 

where /iab(m) = hBA{u) has a non-negative trace, 



S^^hAB{u) > 0, 



(48) 
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but is arbitrary otherwise. For n = 4, any such metric can be interpreted as a 
combined gravitational and electromagnetic plane wave. If equality holds in 
(|48)) . the spacetime is Ricci-flat and, thus, a pure gravitational wave. For any 
choice of /i^s(u), the vector field dy is lightlike and absolutely parallel. For 
a discussion of the geometry of plane-wave spacetimes the reader is refered 
to Penrose (15j . 

Plane-wave spacetimes have the following interesting property, first dis- 
covered by Penrose [16] . Near any lightlike geodesic in any spacetime, the 
metric takes the form of a plane wave in a well-defined limit, called the 
Penrose limit. Thereby the original geodesic is represented as the curve 
{x^ = 0, . . . ,a;"'~^ = 0,v = 0,u) in the limiting plane-wave spacetime (|47)) . 
The Penrose limit can be conveniently written in terms of lightlike Fermi nor- 
mal coordinate, as recently demonstrated by Blau, Frank and Weiss [2]. In 
particular, their analysis showed that Brinkmann coordinates for plane waves 
are lightlike Fermi normal coordinates; in this case all the higher-order terms, 
which are indicated in (j46p by ellipses, vanish exactly. The curvature tensor 
is given along the geodesic {x^ = 0, . . . , a;"~^ = 0,v = 0,u) by 

RAnBn{u) = hAB{u) (49) 



and R^i,vcrT{u) = for all other index combinations, compare p6|l with (j47l) . 

We want to write the generalized Jacobi equation near the lightlike geodesic 
(a;^ — 0, . . . , x"^^ = 0, u = 0, u). We will use u — x"^ for the curve parame- 
ter. After calculating the contravariant metric components from (I47p . we can 
write the Hamiltonian ([23|l : 

H{x,p) ^ Pu ~ \hAB{u)x'^ x^ p^ - - — J^AJP^ _ (50) 
2 2p^ 

Note that this Hamiltonian is of second order with respect to the coordi- 
nates x^ and independent of the coordinate v. Hamilton's equations with 
the Hamiltonian (l50l) take the form 



du _ ^ dv _ 1^ ^A^B S'^^PAPB dx^ _ 5^^ PB 
du ^ du 2 2p1 ^ du Pv ' 

(51) 

dpu_ldhAB AB ^-hA^x^'n (^2) 

— — — - X X Py , — — U, — — UAB X Py . l^OZj 

au I du du du 

Applying d/du to ([^ . and using ([5^ . gives us the lightlike geodesic equation 
in second-order form: 

^--S^^'h.cx^, (53) 
d^v 1 dhAB A B ou A dx^ 



If we write 



du^ - 2 du ^"-^^-^^ABX^^. (54) 
x^{u) = + C^(u) , v{u) = + C{u) , (55) 
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V — u 



Fig. 2 Past light cone of an event in a plane-wave spacetime (|47p . Three dimensions 
{x^,v, u) are shown, with u future-pointing lightlike (from lower right to upper left), 
V future-pointing lightlike (from lower left to upper right) , and orthogonal to 
both. The picture is valid for the case that is an eigendirection of hAB{u) with 
positive eigenvalue. In this case in the three-dimensional picture all the light rays 
issuing from an event into the past are refocused into another event, with the 
exception of one single light ray that stays on an integral curve of the absolutely 
pa ralle l lightlike vector field 9„ . A colour version of this picture ca n b e found online 
in [13. For a similar picture, hand-drawn by Roger Penrose, see [15] . 



and linearize with respect to and we get the generalized Jacobi equa- 
tion 

For the transverse coordinates = f^, these equations are independent 
of the velocities; hence, the generalized Jacobi equation coincides with the 
ordinary geodesic equation. Moreover, the (generalized) Jacobi equation even 
coincides with the exact geodesic equation (|53|) . If we solve these equations 
with initial conditions 

a;'^(wo)=0, -^^B = e^ (57) 

it gives us the cross-section of an initially circular light bundle around the 
geodesic {x^ = 0, ...,a;"~^ = 0,?; = 0, u) with vertex at m = wq. Such a 
bundle will have an elliptic cross-section, for arbitrarily large opening angle. 
The rate of expansion is positive in eigen-directions of Kab (u) with negative 
eigenvalues and negative in eigen-directions of Hab (u) with positive eigenval- 
ues. Condition (|48p makes sure that at least one eigenvalue is positive (unless 
all are zero in which case (|47p is just the Minkowski metric). In directions 
with positive eigenvalue the light bundle will be refocussed until a conjugate 
point is reached. This focusing property is illustrated by Figure [2] which dis- 
plays the past-light cone of an event in a plane-wave spacetime, with three 
dimensions {x^,v,u) shown. 
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The fact that, for the hghthke geodesic {x^ — 0, . . . , x"^^ — 0,v — 0,u) 
in a plane-wave spacetime, the generahzed Jacobi equation coincides with 
the ordinary Jacobi equation has the foUowing interesting consequence: In 
the Penrose hmit the difference between the generahzed Jacobi equation and 
the Jacobi equation vanishes. 



5 Concluding remarks 

In this article we have discussed a generalized Jacobi equation for lightlike 
geodesies, following as closely as possible the approach that was brought 
forward by Hodgkinson, Mashhoon and others for timelike geodesies. As an 
alternative, one could also apply the approach of Bazahski [1], which was 
mentioned in the introduction, to lightlike geodesies. In that case one would 
consider hghthke geodesies of the form x{s) = X{s) + e^(s), where X{s) 
is a lightlike reference geodesic, and solve the geodesic equation for x{s) 
iteratively up to some order N with respect to e. For any finite order TV, the 
resulting equation would not be valid for light bundles of arbitrarily large 
opening angle, in contrast to the generalized Jacobi equation treated here. 
On the other hand, for e sufficiently small it would be valid for arbitrarily 
large parameter intervals. An interesting application of this Bazahski-type 
approach to lightlike geodesies, which apparently has not been considered in 
the literature so far, could be to study the caustics of light bundles up to 
some order N with respect to e. 
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